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Using Monte Carlo simulations on different system sizes we determine with high precision the 
critical thresholds of two families of directed percolation models on a square lattice. The thresholds 
decrease exponentially with the degree of connectivity. We conjecture that p c decays exactly as the 
inverse of the coodination number. 
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Directed percolation (DP) describes generically the dy- 
namics of adsorbing processes and has been applied to 
epidemics, forest fires, surface catalysis, etc [1-6] DP dis- 
plays a phase transition as function of the propagation 
probability p between an "adsorbing" state and a per- 
colating (or active) state at a critical threshold p c . The 
scaling properties of DP have been known since over two 
decades. 

The value of p c depends on the lattice and the rule of 
connectivity. In 1+1 dimension the most studied case is 
the tilted square lattice with nearest neighbor connectiv- 
ity (see Fig. 1) giving p c = 0.6447 ± 0.0001 0. For many 
applications it is of interest to consider a longer range of 
connectivities. In particular also an infinite range model 
with a probability decaying as a powerlaw with distance 
has been studied @. However, what is astonishingly 
missing in the literature are studies of propagation prob- 
abilities of intermediate range although these are in prac- 
tice the most common situations. It is therefore the aim 
of the present paper to calculate the percolation thresh- 
olds of two families of models in 1+1 dimensions having 
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FIG. 1: Schematic plot of the definitions of (a) model I and 
(b) model II. The full lines correspond to (a) k = and (b) 
k — 1, adding the dashed lines one obtains (a) k = 4 and (b) 
k = 3 and adding the dotted lines one obtains (a) k = 6 and 
(b) k = 5. 
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FIG. 2: Double-logarithmic plot of the percolation thresholds 
as function of k for model I (triangles) and of k — 1 for model 
II (circles). 



finite, varying range of interactions. 

Model I is defined on a tilted square lattice as shown 
in Fig. la. The connectivity k is given as the number 
of neighbors to which a site i can be connected in the 
row below. In model I the k is always even. Fig. la also 
shows the cases k — 2,4 and 6. Model II is defined on 
the standard square lattice and in Fig. lb we can see the 
cases k — 1, 3 and 5. For this case, k is always odd. 

We studied the two models up to k = 15 on lattices 
of sizes ranging from 256 x 20000 till 2048 x 20000. We 
simulated for each value of k and L the model for differ- 
ent values of p averaging over 2000 initial configurations 
and monitored if the system percolated or not. From 
the inflection point of the histogram as function of p we 
determined p c {L). Finally we extrapolated to the ther- 
modynamic limit p c through the expression 



Pc (L) =p c (\+aL- x / v ) 



(1) 
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In table 1 we show the obtained thresholds with their 
estimated error bars and in Fig. 2 we see these values in 
a double-logarithmic plot. Since the slope is —1 ± 0.05 
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we conclude that 

p c = r]/k (2a) 
with 77 = 1.3 ± 0.1 for model I and 

Pc = v /(k - 1) (26) 

with T7 = 0.9 ± 0.1 for model II. 

In summary, we have described in this short paper 
numerical evidence that the percolation threshold of di- 
rected percolation decays like rj/k as function of the coor- 
dination k. This result is consistent with the dependence 
of the percolation threshold on the coordination number 
on the Cayley tree |j| and we conjecture it to be exact. 



As a future work, it would be interesting to calculate the 
thresholds for higher dimensional lattices. 
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II 


2 


0.6447 ± 0.0002 


3 


0.4395 ± 0.0003 


4 


0.3272 ± 0.0002 


5 


0.2249 ± 0.0003 


6 


0.2121 ±0.0003 


7 


0.1470 ±0.0002 


8 


0.1553 ±0.0003 


9 


0.1081 ±0.0002 


10 


0.1220 ±0.0002 


11 


0.0851 ±0.0002 


12 


0.0999 ± 0.0002 


13 


0.0701 ± 0.0002 


14 


0.0846 ± 0.0003 


15 


0.0549 ± 0.0002 



TABLE I: Percolation threshold of directed percolation to 
model I and model II. 
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